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Abstract 

We consider power allocation algorithms for fixed-rate transmission over Nakagami-m non- 
ergodic block-fading channels with perfect transmitter and receiver channel state information and 
discrete input signal constellations, under both short- and long-term power constraints. Optimal 
power allocation schemes are shown to be direct applications of previous results in the literature. 
We show that the SNR exponent of the optimal short-term scheme is given by m times the Singleton 
bound. We also illustrate the significant gains available by employing long-term power constraints. 
In particular, we analyze the optimal long-term solution, showing that zero outage can be achieved 
provided that the corresponding short-term SNR exponent with the same system parameters is strictly 
greater than one. Conversely, if the short-term SNR exponent is smaller than one, we show that 
zero outage cannot be achieved. In this case, we derive the corresponding long-term SNR exponent 
as a function of the Singleton bound. Due to the nature of the expressions involved, the complexity 
of optimal schemes may be prohibitive for system implementation. We therefore propose simple 
sub-optimal power allocation schemes whose outage probability performance is very close to the 
minimum outage probability obtained by optimal schemes. We also show the applicability of these 
techniques to practical systems employing orthogonal frequency division multiplexing. 
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I. Introduction 

A key design challenge for wireless communications systems is to provide high-data-rate 
wireless access, while optimizing the use of limited resources such as available frequency 
bandwidth, transmission power and computational ability of portable devices. Reliable trans- 
mission is particular challenging for wireless communications systems due to the harsh, 
time- varying signal propagation environment. Mobility and multipath propagation [1], [2], 
[3] lead to time- selective and frequency selective fading channels, where the dynamics of the 
signal variations depend on mobile velocity, carrier frequency, transmission bandwidth, and 
the particular scattering environment. 

The use of orthogonal frequency division multiplexing (OFDM) technologies is a proven 
approach for providing high data rates in wireless communications systems. Standards such 
as IEEE 802.11 (WiFi) [4] and IEEE 802.16 (WiMax) [5] already include OFDM as a core 
technology, and future generations of mobile cellular systems are likely to also feature multi- 
carrier techniques. OFDM transmission over frequency- selective or time-frequency-selective 
wireless fading channels is adequately modelled as a block-fading channel. 

The block-fading channel [6], [1] is a useful channel model for a class of time- and/or 
frequency-varying fading channels where the duration of a block-fading period is determined 
by the product of the channel coherence bandwidth and the channel coherence time [7]. Within 
a block-fading period, the channel fading gain remains constant, while between periods the 
channel gains change according to a system-specific rule. In this setting, transmission typically 
extends over multiple block-fading periods. Frequency-hopping schemes as encountered in the 
Global System for Mobile Communication (GSM) and the Enhanced Data GSM Environment 
(EDGE), as well as transmission schemes based on multiple antenna systems, can also 
conveniently be modelled as block-fading channels. The simplified model is mathematically 
tractable, while still capturing the essential features of the practical transmission schemes 
over fading channels. 

In many situations of practical interest, channel state information (CSI), namely the de- 
gree of knowledge that either the transmitter, the receiver, or both, have about the channel 
gains, greatly influences system design and performance. In general, optimal transmission 
strategies over a block-fading channel depend on the availability of CSI at both sides of 
the transmission link [1]. At the receiver side, time- varying channel parameters can often be 
accurately estimated [7]. Thus, perfect CSI at the receiver (CSIR) is a common and reasonable 
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assumption. Conversely, perfect CSI at the transmitter (CSIT) depends on the specific system 
architecture. In a system with time-division duplex (TDD), the same channel can be used for 
both transmission and reception, provided that the channel varies slowly. In this case, perfect 
CSIR can be used reciprocally as perfect CSIT [8]. In other system architectures, CSIT 
is provided through channel- state-feedback from the receiver. When no CSIT is available, 
transmit power is commonly allocated uniformly over the blocks. In contrast, when CSIT 
is available, the transmitter can adapt the transmission mode (transmission power, data rate, 
modulation and coding) to the instantaneous channel characteristics, leading to significant 
performance improvements. 

We distinguish between two cases of transmission dynamics. On the one hand, if no delay 
constraints are enforced, transmission extends over a large (infinite) number of fading blocks. 
The corresponding fading process is stationary and ergodic, revealing the fading statistics 
during the transmission. The maximum data rate for this case, termed the ergodic capacity, 
was determined in [9], assuming perfect CSI at both transmitter and receiver. Two coding 
schemes have been shown to achieve the ergodic capacity. In [9], a variable-rate, variable- 
power transmission strategy based on a library of codebooks, and driven by the CSIT, was 
suggested. In contrast, a fixed-rate, variable-power transmission strategy was proposed in [10] 
based on a single codebook, providing a practically more appealing alternative in the form 
of a conventional Gaussian encoder followed by power allocation driven by the CSIT. 

On the other hand, when a delay constraint is enforced, the transmission of a codeword 
only spans a finite number of fading blocks. This constraint corresponds to real-time transmis- 
sion over slowly varying channels. Therefore, this situation is relevant for wireless OFDM 
applications in wireless local area networks (WLAN). As the channel relies on particular 
realizations of the finite number of independent fading coefficients, the channel is non-ergodic 
and therefore not information stable [11], [12]. It follows that the Shannon capacity under 
most common fading statistics is zero, since there is an irreducible probability, denoted 
as the outage probability, that the channel is unable to support the actual data rate, [6], 
[1]. For sufficiently long codes, the word error rate is strictly lower-bounded by the outage 
probability. In some cases there is a maximum non-zero rate and a minimum finite signal- 
to-noise ratio (SNR) for which the minimum outage probability is zero. This maximum rate 
is commonly referred to as the delay-limited capacity [13]. In this paper, we will consider 
fixed-rate transmission strategies over delay-limited non-ergodic block-fading channels. 

In a practical system, only causal CSIT is available. Thus, in general, the channel gains 
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are only known up to (and possibly including) the current block-fading period. However, 
in an OFDM system with multiple parallel carriers, the causal constraint still allows for all 
sub-carrier channel gains to be known simultaneously in a seemingly non-causal manner, as 
compared to a block-fading channel based on frequency-hopping single-carrier transmission. 
Here, we will only consider the OFDM-inspired scenario where perfect non-causal CSIT is 
available. 

As mentioned above, when perfect CSI is available at the transmitter, power allocation 
techniques can be used to increase the instantaneous mutual information, thus improving the 
outage performance. Multiple power allocation rules derived under a variety of constraints 
have been proposed in the literature [12], [14], [15], [16], [17]. The optimal power allocation 
minimizes the outage probability subject to a short-term power constraint over a single 
codeword or a long-term power constraint over all transmitted codewords. The optimal 
transmission strategy, subject to a short-term power constraint, was shown in [12] to consist of 
a random code with independent, identically distributed Gaussian code symbols, followed by 
optimal power allocation based on water- filling [18]. The optimal power allocation problem 
is also solved in [12] under a long-term power constraint, showing that remarkable gains are 
possible with respect to transmission schemes with short-term power constraints. In some 
cases, the optimal power-allocation scheme can even eliminate outages, leading to a minimum 
outage probability approaching zero [12], [19], and thus, a non-zero delay-limited capacity. 
In particular, gains of more than 12 dB are possible at practically relevant error probabilities. 
Again, the optimal input distribution is Gaussian. The optimal power allocation problem 
under a long-term power constraint, and with perfect CSIR but only partial CSI available at 
the transmitter is considered in [20]. The problem is solved for the limiting case of large 
SNR, leading to similar impressive improvements in outage performance. 

In practical wireless communications systems, coding schemes are constructed over discrete 
signal constellations, e.g., PSK, QAM. It is therefore of practical interest to derive power 
allocation rules for coded modulation schemes with discrete input constellations, minimizing 
the outage probability. A significant step towards this goal was achieved in [21], where the 
fundamental relationship between mutual information and MMSE developed in [22] proved 
instrumental to optimizing the transmit power of parallel channels with discrete inputs. As 
stated in [21], the developed power allocation rule for parallel channels can be applied directly 
to minimize the outage probability of delay-limited block-fading channels under short-term 
power constraints. However, the optimal solution in [21] does not reveal the impact of the 
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system parameters involved, and may also be prohibitively complex for practical applications 
with computational power and memory limitations. 

In this paper, we study power allocation schemes that minimize the outage probability of 
fixed-rate coded modulation schemes using discrete signal constellations under short- and 
long-term power constraints. In particular, we study classical coded modulation schemes, as 
well as bit-interleaved coded modulation (BICM) using suboptimal non-iterative decoding 
[23]. Similarly to the uniform power allocation case, we show that under a short-term power 
constraint, an application of the Singleton bound [24], [25], [26], [27] leads to the optimal 
SNR exponent (diversity gain) of the channel. In particular, we show that for Nakagami-m 
channels, the optimal SNR exponent is given by m times the Singleton bound [24], [25], 
[26], [27]. In the long-term case, we derive the optimal power allocation scheme. We show 
that the underlying structure of the solution for Gaussian inputs in [12] remains valid, where 
no power is allocated to bad fading realizations, minimizing power wastage. We also show 
that the relationship between the mutual information and the minimum- mean- squared-error 
(MMSE) reported in [22] is instrumental in deriving the optimal outage-minimizing long- 
term solution^ We analyze the optimal long-term solution, showing that zero outage can be 
achieved provided that the SNR exponent corresponding to the short-term scheme with the 
same parameters is strictly greater than one, implying the delay-limited capacity is non-zero. 
Conversely, if the short-term SNR exponent is smaller than one, we show that zero outage 
cannot be achieved. In this case, we derive the corresponding long-term SNR exponent as a 
function of the Singleton bound. 

Practical transmitters may have limited memory and computational resources that may 
prevent the use of the optimal solution based on the MMSE. We further aim at reducing 
the computational complexity and memory requirements of optimal schemes by proposing 
sub-optimal short- and long-term power allocation schemes. The sub-optimal schemes enjoy 
significant reductions in complexity, yet they only suffer marginal performance losses as 
compared to relevant optimal schemes. For the suggested sub-optimal schemes, we further 
characterize the corresponding SNR exponents and delay-limited capacities for short- and 
long-term constraints as functions of the modified Singleton bound. 

The paper is further organized as follows. In Section [IlJ the system model, basic as- 

'The optimal power allocation algorithm with discrete inputs has been independently reported in [28]. We became aware 
of the results in [28] after the submission of the conference version of this paper [29]. 



February 5, 2008 



DRAFT 



6 



sumptions and related notation are described, while mutual information, MMSE and outage 
probability are introduced in Section [jllj Optimal and sub-optimal power allocation schemes 
with short-term power constraints are considered in Section |IVJ and corresponding optimal 
and sub-optimal power allocation schemes with long-term power constraints are investigated 
in Section |Vj Numerical examples are used throughout the paper to illustrate the presented 
results. Concluding remarks, summarized in Section [VTJ complete the main body of the paper. 
To support the readability of the paper, lengthy proofs are moved to appendices. 

Throughout the paper, we shall make use of the following notation. Scalar and vector 
variables are characterized with lowercase and boldfaced lowercase letters, respectively. 
The expectation with respect to the fading statistics is simply denoted by E [■], while the 
expectation with respect to any other arbitrary random variable $ is denoted by E$ [•]. Fur- 
thermore, the expectation with respect to an arbitrary random variable $, with the constraint 
$ E TZ is denoted as E^ en [■]. We define (x) = -| Ya= 

t X{ clS the arithmetic mean of x = 
(xi, X2, xb)- The exponential equality /(£) = K^~ d indicates that lim^oo — K, 

with the exponential inequalities <,> similarly defined. IR™ = {( G l n |( > 0}, min{a, b} 
denotes the minimum of a and b, \£] (|_£J) denotes the smallest (largest) integer greater 
(smaller) than f, while {f{x)) + = if f(x) < 0, and {f{x)) + = f(x) if f(x) > 0. 

II. System Model 

Consider transmission over an additive white Gaussian noise (AWGN) block-fading channel 
with B blocks of L channel uses each, in which, for b — 1, . . . , B, block b is affected by a 
flat fading coefficient /ij G C. The corresponding block diagram is shown in Figure [T] Let 
lb = \hb\ 2 be the power fading gain and assume that the fading gain vector 7 = (71, . . . , 7b) 
is available at both the transmitter and the receiver. The transmit power is allocated to the 
blocks according to the scheme p(7) = (^1(7), • • • ,Pb(t))- Then, the complex baseband 
channel model can be written as 

yb = Vpb~h)h b x b + z b , b=l,...,B, (1) 

where y b E C L is the received signal in block b, x b E X L C C L is the portion of the 
codeword being transmitted in block b, X C C is the signal constellation and z b E C L is a 
noise vector with independent, identically distributed (i.i.d.) circularly symmetric Gaussian 
entries ~ A/"c(0, 1). Assume that the signal constellation X is normalized in energy such that 
J2x&x \ x \ 2 = 2 M > where M = log 2 \ X\. Then, the instantaneous received SNR at block b is 
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given by Pb{~{)lb- The following power constraints are considered [12] 



1 B 

Short-Term: (p( 7 )> 4 _ ^p 6 ( 7 ) < P (2) 

6=1 



Long-Term: E [(p( 7 ))] = E 



1 B 



B 

b=l 



< P. (3) 



In all cases, P represents the average SNR at the receiver. We will denote by p st and 
Pit the power allocation vectors corresponding to short- and long-term power constraints, 
respectively. We will also denote by p eq (p) = (p, ■ ■■,p) the uniform power vector that 
allocates the same power p to each block. 

We consider block-fading channels where hb are independent realizations of a random 
variable H, whose magnitude is Nakagami-m distributed [30], [31] and has a uniformly 
distributed phase^J The fading magnitude has the following probability density function (pdf) 



2m m h 



ml,2m-l 



f\ H \(h) = ""I" e~ mh \ (4) 
1 [m) 

where T(a) is the Gamma function, T(a) = t a ~ 1 e~ t dt [32]. The coefficients % are 
realizations of the random variable \H\ 2 whose pdf and cdf are given by 

m m m— 1 

Ul) = < " (5) 




and 

P, b (0={ L(m > ~ (6) 

respectively, where T(a, £) is the upper incomplete Gamma function, T(a,£) = t a ~ 1 e~ t dt 
[32]. 

The Nakagami-m distribution encompasses many fading distributions of interest. In par- 
ticular, we obtain Rayleigh fading by letting m = 1 and an accurate approximation of Rician 
fading with parameter K by setting m = (K + 1) 2 /(2K + 1) [31]. 

2 Due to our perfect transmitter and receiver CSI assumption, we can assume that the phase has been perfectly compensated 
for. 



February 5, 2008 



DRAFT 



s 



III. Mutual Information, MMSE and Outage Probability 

The channel model described in ([T]) corresponds to a parallel channel model, where 
each sub-channel is used a fraction ^ of the total number of channel uses per codeword. 
Therefore, for any given power fading gain realization 7 and power allocation scheme p(7), 
the instantaneous input-output mutual information of the channel is given by [18] 

1 B 

Jb(p(7)> 7) = ^ Yl J *(^ fe )' ( 7 ) 
6=1 

where Ix(p) is the input-output mutual information of an AWGN channel with input con- 
stellation X and received SNR p. In this paper we will consider coded modulation (CM) 
schemes with uniform input distribution, for which Ix{p) is given by 



jCM 



L X 



(P) 



ll) g 2 ( Yl e 

Furthermore, we also consider bit-interleaved coded modulation (BICM) using the classical 
sub-optimal non-iterative BICM decoder proposed in [33], for which the mutual information 
for a given binary labeling rule^jcan be expressed as the mutual information of M binary-input 
continuous-output-symmetric parallel channels [23], 

e -\Jp{X-x')+Z\> 



-\^p(x-x')+Zf+\Z\ 2 



1 M 



/-BICM 
l X 



(p) 



9=0 j=l x& xl 



log 2 



E 

x'ex 



E 

x'dX 3 a 



-\^p(X-x')+Z\ 2 



where the sets X l c contain all signal constellation points with bit c in the j-th binary labeling 
position. Both CM and BICM mutual information expressions can be efficiently evaluated 
numerically using Gauss-Hermite quadratures [32]. 

A fundamental relationship between the MMSE and the mutual information (in bits) in 
additive Gaussian channels is introduced in [22] showing that 



^Ix(p) = -^-MMSE^(p), 
dp log 2 



(8) 



where MMSE^-(p) is the MMSE for a given signal constellation X expressed as a function 
of the SNR p. This relationship proves instrumental in obtaining optimal power control rules. 



3 We select Gray labeling [34], since it has been shown to maximize the mutual information for the non-iterative BICM 
decoder [23]. 
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In particular, for the CM case, the MMSE resulting from estimating the input based on the 
received signal over an AWGN channel with SNR p can be written as [21], 

MMSE£ M (p) = E \\X - X\ 2 ~\ = E [\X - E [X\Y] | 2 ] (9) 



i- 1 



2 



-dy, (10) 



where X = E [X\Y] is the MMSE estimate of X given the channel observation Y. Once 



again, ( |I0| ) can be efficiently evaluated using Gauss-Hermite quadratures [32]. In the case 
of BICM, we obtain an equivalent set of symmetric binary-input continuous output channels 
[23]. However, due to the demodulation process, the equivalent channels have a noise that is 
non-Gaussian, and more importantly, non-additive. We therefore define the function derivative 
of lf CM (p), denoted by MMSE| ICM (p), by enforcing the relationship © to hold as follows 

MMSE| ICM (p) 4 ^lf CM {p). (11) 
dp 

Note that this is only a shorthand notation, so that whenever MMSE^ appears in the coming 
sections, it can be replaced by either MMSE£ M (p) or MMSE^, ICM (p). Therefore, ([n} does 



not denote the MMSE in estimating the input bits given the noisy channel observation. The 
function MMSE| ICM (p) can again be easily evaluated numerically. 

Finally, we define the transmission to be in outage when the instantaneous input-output 
mutual information is less than the target fixed transmission rate R. For a given power 
allocation scheme p(7) with power constraint P, the outage probability at transmission rate 
R is given by [6], [1] 

P out (p( 7 ),P,R) = Pr(/ B (p( 7 ),7) < R) 



Pr (^E J *(^) <i? )- w 



Since J^ ICM (p) < I^ M (p) we will have that the corresponding outage probabilities verify 
that P™ M (p(j),P,R) > P^(p(j),P,R). All the algorithms and results presented in the 
following are valid for both CM and BICM. Therefore, unless explicitly stated, we will use 
the common notation Ix{p) and MMSE^(p) to refer to both. 

IV. Short-Term Power Allocation 

Short-term power allocation schemes are applied to systems where the transmit power of 
each codeword is limited to BP. Following the definition of short-term power constraint 
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given in Section[TTJ a given short-term power allocation scheme p st = (p%, . . . ,Pb) must then 
satisfy | J2f=i Vb < P. 

A. Optimal Short-Term Power Allocation 

The optimal short-term power allocation rule p° pt (7) is the solution to the outage proba- 
bility minimization problem [12]. Mathematically we express p° pt (7) as 

P°t Pt (7) = arg min P out (p,P,R). (13) 

For short-term power allocation, the power allocation scheme that maximizes the instanta- 
neous mutual information at each channel realization also minimizes the outage probability 
since the available power can only be distributed within one codeword. Formally, we have 
[12] 

Lemma 1: Let p° pt (7) be a solution of the problem 

Maximize Ef=i -MP&7&) 
< Subject to iJ2tiPb<P < 14 ) 
p b >0,b=l,...,B. 

Then p° t pt (7) is a solution of flD) . 

Proof: See Appendix [j] ■ 
The solution of problem ( [14] ), which is based on the relationship between the MMSE and 
the mutual information [22], was obtained in [21]. From [21] one has the following theorem. 



Theorem 1: The solution of problem ([14]) is given by 



p°pt( 7 ) = i-MMSE^ 1 Lin jl, l\\ , (15) 
for b — 1, . . . , B, where rj is chosen such that the power constraint is satisfied, 

itlMMSE-( m i„{l,^})=P. (16, 

Proof: See [21] for details. ■ 
As outlined in [21], the results of Theorem [TJ are valid for any constellation. In particular, 
since the MMSE for Gaussian inputs (GI) is 

MMSE °' W = icg2(TTrt < <17) 
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the inverse function can be written in closed form [21] and we therefore recover the waterfillng 
solution by replacing MMSE^(p) by MMSE^(p) [12]. 

The optimal short-term power allocation scheme improves the outage performance of coded 
modulation schemes over block-fading channels. However, it does not increase the outage 
diversity compared to a uniform power allocation, as shown in the following result. 

Proposition 1: Consider transmission over the block-fading channel defined in ([T]) with 
the optimal power allocation scheme p° t pt (7) given in ([15]). Assume input constellation size 
\X\ = 2 M . Further assume that the power fading gains follow the distribution given in ([5]). 
Then, for large P and some /C opt > the outage probability behaves as 

PoMfh), p, R) = ic opt p- mdBiR) , (is) 

where d B (R) is the Singleton bound given by 

R 



d B (R) = 1 + 
Proof: See Appendix [IJ 



B 1 

M 



(19) 



B. Sub-optimal Short-Term Power Allocation Schemes 

Although the power allocation scheme in ( [l"5] ) is optimal, it involves an inverse MMSE 
function, which may be excessively complex to implement or store for specific low-cost 
systems. Moreover, the MMSE function provides little insight into the role of each system 
parameter. In this section, we propose sub-optimal power allocation schemes similar to water- 
filling that tackle both drawbacks, leading to only minor losses in outage performance as 
compared to the optimal solution. 



1) Truncated water-filling scheme: The complexity of the solution in ( fl5| ) is due to the 
complex expression of Ix(p) m problem ( [T4| ). Therefore, in order to obtain a simple sub- 
optimal solution, we propose an approximation for Ix{p) in problem ([14]). For Gaussian 
input channels with J(p) = log 2 (l + p), optimal power allocation is obtained by the simple 
water- filling scheme [18]. This suggests the use of the following approximation for Ix(p)- 

flog,(l + p), p<p 
(20) 
log 2 (l + /3), otherwise, 



February 5, 2008 



DRAFT 



12 



where (3 is a design parameter to be optimized for best performance. An example of I tw (p) 
is illustrated in Figure [2] The resulting sub-optimal scheme p^i'j) is given as a solution of 



Maximize J2b=i Jtw 



{Pblb) 



Subject to ^T,b=iPb< p 

Pb>0, b=l,...,B. 



(21) 



Theorem 2: A solution to the problem in ( |2T| ) is given by 

11 _B 2^6=1 7t , - ^ 



K w (t) 



it. ' 



mm 



V 



otherwise, 



(22) 



for b = 1, . . . , B, where r\ is chosen such that 

B 

mil) < — . | r] 



1 S 

— > mil 
B ^ 

6=1 



I 
lb' 



1 

7fc 



P. 



(23) 



Proof: See Appendix |Ilj 
Without loss of generality, assume that 71 > . . . > 7s, then, similarly to water-filling, r] 
can be determined such that [12] 



(k-i) V = BP-y2 — + E 1 



(24) 



where k, I are integers satisfying — < i] < 



6=1 6— J 

and ^ 



< ?y < 



From Theorem |2} the resulting power allocation scheme is similar to water-filling, except 
for the truncation of the allocated power at — . We refer to this scheme as truncated water- 
filling. The outage performance obtained by the truncated water-filling scheme depends on 
the choice of the design parameter (3. We now analyze the asymptotic performance of the 
outage probability, thus providing some guidance for the choice of [3. 

Proposition 2: Consider transmission over the block-fading channel defined in ([T]) with 
input constellation X and the truncated water- filling power allocation scheme P^t) given 
in ( |22| ). Assume that the power fading gains follow the distribution given in ([5]). Then, for 
large P, the outage probability P out (p*"(7), P, R) is asymptotically upper bounded by 



Pouter (7),^) <}CpP- md ^ R \ 



where 



dfs(R) = 1 + 



B 1 



R 



(25) 



(26) 



and Ix((3) is the input-output mutual information of an AWGN channel with SNR (3. 
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Proof: See Appendix |nij ■ 
From the results of Proposition [TJ Proposition [5J and noting that P out (pg™(7), P, R) > 
^out(Pst Pt (7),^^)> we have 

P out (p^(7), P, R) = JC tw P~ md ^ R \ (27) 

where d tw (R) satisfies that dp(R) < d tw (R) < d B (R). Therefore, the truncated water- filling 
scheme is guaranteed to obtain optimal diversity whenever dp(R) = d B (R), or equivalently, 
when 

R 



Bll -M 



(28) 



which implies that 



13 - '* I B-[B(l-f,)\ 1 ^ & 



Therefore, by letting (3 — > oo, the truncated water-filling power allocation scheme given in 
( |2~2~1 ), which now becomes the classical water-filling algorithm for Gaussian inputs, provides 
optimal outage diversity at any transmission rate. For any rate R that is not at a discontinuity 
point of the Singleton bound, i.e. R such that B (l — ||j is not an integer, we can always 
design a truncated water-filling scheme that obtains optimal diversity by choosing (3 > (3r. 

With the results above, we choose (3 as follows. For a transmission rate R that is not a 
discontinuity point of the Singleton bound, we perform a simulation to compute the outage 
probability obtained by truncated water-filling with various [3 > (3r and pick the (3 that gives 
the best outage performance. The dashed lines in Figure [3] illustrate the performance of the 
obtained schemes for block-fading channels with B = 4 and QPSK input under Rayleigh 
fading. At all rates of interest, the truncated water-filling schemes suffer only minor losses 
in outage performance as compared to the optimal schemes (solid lines), especially at high 
SNR. We also observe a remarkable difference with respect to pure water-filling for Gaussian 
inputs (dotted lines). As a matter of fact, pure water-filling performs worse than uniform 
power allocation. 

For rates at the discontinuous points of the Singleton bound, especially when operating 
at high SNR, [3 needs to be relatively large in order to maintain diversity. However, large (3 
increases the gap between I tw (p) and Ix(p), thus degrading the performance of the truncated 
water-filling scheme. For (3 = 15, the gap is illustrated by the dashed lines in Figure |4j In the 
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extreme case where (3 — > oo, the truncated water-filling turns into the water-filling scheme, 
which exhibits a significant loss in outage performance as illustrated in Figure [3] 

2) Refined truncated water-filling schemes: To obtain a better approximation to the optimal 



power allocation scheme, we need a more accurate approximation to Ix(p) in ( |T4| ). We 
propose the following bound. 



i x { P )<r et {p) 



log 2 (l + p), P<« 

fdog 2 (p) + a, a < p < (3 (30) 
/clog 2 (/3) + a, otherwise, 
where « and a are chosen such that (in a dB scale) /tlog 2 (p) + a is a tangent to Ix(p) at 
a predetermined point p . Therefore a is chosen such that n\og 2 (a) + a = log 2 (l + a), 
and (3 is a design parameter. These parameters are reported in Table [j] for CM and BICM 
using various modulation schemes. An example of the approximation is also illustrated by 
the dashed-dotted curve in Figure [2} 

The corresponding optimization problem can be written as 



Maximize YjLi ^(P&Tfc) 



Subject to Ylb=i Pb < BP 

Pb >0,b=l,...,B. 

The refined truncated water- filling scheme Ps£ f (7) is given by the following theorem. 



(3D 



Theorem 3: A solution to problem pT) is 

P 



lb 



1,...,B 



(32) 



if i 2~2b=i ~ < P> anc ^ otherwise, 



Pf 



it " 

KT], 
a 

lb " 

v - 

0. 



> 



«7b 



K7b — ' K 



Q + l 



76 



76 ' 



< 77 < ^ 
lb — ' K lb 

J_ < „ < «±1 
lb — 1 lb 

otherwise, 



(33) 



for b 



B, where r\ is chosen such that 



1 B 
B^2 

6=1 



Pf 



P. 



(34) 
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Proof: See Appendix [TV} 



The refined truncated water-filling scheme provides significant gain over the truncated 
water-filling scheme, especially when the transmission rate requires a relatively large (3 to 



mance of the refined truncated water-filling scheme for block-fading channels with B = 4, 
and QPSK input under Rayleigh fading. The outage performance of the refined truncated 
water-filling scheme is close to the outage performance of the optimal case even at the rates 
where the Singleton bound is discontinuous, i.e. rates R = 0.5, 1.0, 1.5. The performance 
gains of the refined scheme over the truncated water-filling scheme at other rates are also 
illustrated by the dashed-dotted lines in Figure [3j 



We consider systems with long-term power constraints, in which the expectation of the 
power allocated to each block (over infinitely many codewords) does not exceed P. This 
problem has been investigated in [12] for block-fading channels with Gaussian inputs. In this 
section, we obtain similar results for channels with discrete inputs, and propose sub-optimal 
schemes that reduce the complexity of the algorithm. 

A. Optimal Long-Term Power Allocation 

Following [12], the problem can be formulated as 



where (p) = sJ2 b =iPb- 

The following theorem shows that the structure of the optimal long-term solution p° t pt (7) 
of [12] for Gaussian inputs is generalized to the discrete-input case. 

Theorem 4: Consider transmission over the block-fading channel given in ([T]) with input 
constellation X. Assume that the power fading gains in 7 follow the distribution given in 
(|5]). Then, the optimal power allocation scheme for systems with long-term constraint P, is 
given by 



maintain the outage diversity. The dashed-dotted lines in Figure [4] show the outage perfor- 



V. Long-Term Power Allocation 




(35) 




otherwise, 




(36) 



where 
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1) P° pt (l) is the solution of the following optimization problem: 



Minimize (p("y)) 



Subject to -J J2 b=1 Ix(pblb) > R 
p b >0,b=l,...,B. 



(37) 



2) s satisfies 



s = oo 

P =E 7e7 ^ ) [( f p°e t ( 7 )>] 



if lim M E 7e(s) [(^P t ( 7 ))]<P 
otherwise, 



in which, 



(38) 



(39) 



?l(^{ 7 G^:(p opt (7)}< S }. 
Proof: See Appendix |V| ■ 
Theorem 5: Consider transmission at rate R over the block-fading channel given in ([T]) 
with input constellation X and realized power fading gain 7, the power allocation scheme 



that minimizes the input power, which is a solution of the problem given in ( [37] ), is given by 



where i] is chosen such that 



1 

76 



MMSE^ 1 ( min 



1 B ( 

s£ J * ( MMSE ^ 



6=1 



1 



R. 



(40) 



(41) 



Proof: See Appendix [VI| 
As in the Gaussian input case [12], the optimal power allocation scheme either transmits with 
the minimum power that enables transmission at the target rate using an underlying dual short- 
term scheme p opt {"i) with short-term constraint (p opt (7)) < s, or turns off transmission 
(allocating zero power) when the channel realization is bad. Therefore, there is no power 
wastage on outage events. 

The solid lines in Figures [5] and [6] illustrate the outage performance of optimal long-term 
power allocation schemes for transmission over 4-block block-fading channels with QPSK 
and 16-QAM inputs and Rayleigh fading (m = 1). The simulation results suggest that for 
transmission rates where ds(R) > 1, zero outage probability can be obtained with finite 
power. This agrees with the results obtained for block-fading channels with Gaussian inputs 
[12], where only for B > 1 zero outage is possible. 
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To provide more insight into this effect, consider the following long-term power allocation 
scheme, 

. . Jp(7), ien(s) 

Pit (7) = < (42) 
I 0, otherwise, 

where £?( 7 ) is an arbitrary underlying short-term power allocation scheme, 

K(s) = { 7 GMf: (p( 7 )> < s} , (43) 
and s is chosen to satisfy the long-term power constraint, 

E [<p lt (7)>] = E 7e7 e (s) [<p( 7 )>] = P- (44) 

Assuming that for any s, p( 7 ) satisfies -| Y^b=\Ix{pblb) > R f° r an Y 7 ^ 7^(s), then the 
resulting outage probability of the corresponding long-term power allocation is 

Pout(pit(7), P(s), R) = Pr( 7 i K(s)). (45) 



For any long-term power constraint P, the long-term power allocation scheme in ((42J) depends 
on the threshold s defined in ( |44] ). Conversely, for any choice of the threshold s, the long-term 
power P(s) is given by 

P{s) = E Jen(s) [(p( 7 ))] . (46) 

We now consider the behavior of the average power P(s) and the corresponding outage prob- 
ability Pout (pit ( 7 )> P( s )i R) when s — >• 00. In particular, we study the long-term exponent 
defined as 

1 / 7~>\ A r lQgPout(Plt( 7 ),P(g),-R) 

d lt (R) = lim —— . (47) 

P(s)^oo log P(s) 

Firstly, consider the following asymptotic relationship between s and P(s). 

Proposition 3: Consider transmission over a block-fading channel with a long-term power 

allocation scheme corresponding to an arbitrary underlying short-term scheme p( 7 ), a thresh- 



old s given in ((42]), and a long-term power constraint P(s) given in ( [46] ). Assume that £3(7) 
is chosen such that asymptotically in s, the outage probability satisfies 

Pout(Pit( 7 ), P(s), R) = Ks~ d W (48) 
for some finite d(R) > 0. Then, 

^-P(s) = fCd(R)s~ di ^ R) (49) 
as 



Proof: See Appendix |VII 
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From the previous proposition, we obtain the following result. 

Theorem 6: Consider transmission over a block-fading channel with a long-term power al- 
location scheme corresponding to an arbitrary underlying short-term scheme ft? (7), a threshold 



s given in < |42] ), and a long-term power constraint P(s) given in ( |46| ). Assume that ft? (7) is 
chosen such that asymptotically in s, the outage probability satisfies 

^out(pit(7), P(s), R) = lCs-*™ (50) 

for some finite d(R) > 0. Then, if d(R) > 1 we have that 

lim P(s) = P th < 00 and d it (R) = 00, 

s— >oo 

while if d(R) < 1, 

d «W - T%j- Pi) 



Proof: See Appendix |VII 



The previous results highlight the effect of the power constraint on the outage performance 
obtained by a specific power allocation scheme pit (7). In particular, P ut(pit(7), P( s ), R) is 
the outage probability of the block-fading channel with power allocation scheme ft?(7) and 
short-term power constraint (ft?(7)) < s, and d(R) is the corresponding short-term outage 
diversity. When a long-term power constraint is applied, the long-term outage diversity d\ t (R) 
is affected in the following way: 

. If d(R) > 1, then d\ t (R) = 00 and Hindoo P(s) = P t h < 00. Therefore, since 

limP out ( Plt (7),P( S ),i?) = 0, (52) 

s^oo 

there exists a threshold long-term power constraint P th beyond which strictly zero outage 
probability is achieved, proving that vanishing error probability can be achieved and that 
reliable communication is possible at rates below R. R is therefore a lower bound to 
the delay limited capacity [13] of the block-fading channel with power constraint P th . 



• If d(R) < 1, then pi] ) gives the relationship between the long- and short-term outage 
diversity. 

In order to apply the previous theorem to analyze the asymptotic behaviour of the optimal 
power allocation for systems with long-term constraints, we consider the following duality 
between p° t pt (7) and ft? opt (7). 

Proposition 4: Consider transmission at rate R over the block-fading channel given in 
([T]) with the optimal long-term power allocation scheme p° t pt (7) given in < |36] ) and a long- 
term power constraint P(s). Then, independent of the fading statistics, the outage probability 
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satisfies 

Poutfcif (7), P{s), R) = Pout(pTh), s, R), (53) 
where p° pt (7) is the optimal power allocation scheme with short-term power constraint 

<PST(7)> < »• 

Proof: See Appendix |VII ■ 



From Theorem [6} we have the following result. 

Corollary 1: Consider transmission at rate R over the block-fading channel given in ([T]) 
with the optimal long-term power allocation scheme p° pt (7). Assume input constellation 
X of size \X\ = 2 M . Further assume that the power fading gain 7 follows a Nakagami-m 
distribution given in (|5]). Then, the delay-limited capacity is non-zero whenever cIb(R) > —■ 
Conversely, when cIb(R) < — , the outage probability asymptotically behaves as 

PL(j%*(rr),P{s), R) = icTP~ dTiR \ (54) 

where P is the long-term power constraint, and <i t pt (i?) is the optimal long-term outage 



diversity given by 



Proof: From Propositions [T] and [4j we have 

PL(pT (7), P{*), R) = lCo P tS- md »M. (56) 

Therefore, the corollary can be obtained as a direct application of Theorem [6} ■ 
This behavior is illustrated in Figure [7J where the outage probability with QPSK inputs, 
m = 0.5 and R — 1.7 has been plotted as a function of the average long-term power P(s) 
and as a function of the dual short-term constraint s. As predicted by the previous results, 
the dual short-term curve has slope mdB(R) =0.5. Furthermore, since ds(R) = 1 < — = 2, 
we observe that the long-term outage curve (as a function of P(s)) has slope <i° t pt (i?) = 

md B (R) _ -y 



l-md B (R) 

B. Sub-optimal Long-Term Power Allocation 

In the optimal long-term power allocation scheme p° t pt (7) given in Theorem |4j s can be 
evaluated offline for any fading distribution. Therefore, given an allocation scheme £? opt (7), 
the complexity required to evaluate p° t pt (7) is low. Thus, the complexity of the long-term 
power allocation scheme is mainly due to the complexity of evaluating p opt (7), which 
requires the evaluation or storage of MMSE^(p) and Ix(p). In this section, we propose 
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sub-optimal long-term power allocation schemes by replacing the optimal underlying short- 
term algorithm p opt (7) with simpler allocation rules. 

A long-term power allocation scheme pit (7) corresponding to an arbitrary p(i) is obtained 
by replacing £> opt (7) in ([36]), ( |3~8] ) and ([39]) with p(j). From ( [36] ), the long-term power 
allocation scheme pi t (7) satisfies 

E [<pi t (7)>] =^yen(s) [(p(7))] = P- (57) 

Therefore, a long-term power allocation scheme corresponding to an arbitrary £3(7) is sub- 
optimal with respect to p° t pt (7). Following the transmission strategy in the optimal scheme, 
we consider the power allocation schemes p(j) that satisfy the rate constraint I B (p{p(), 7) > 
R to avoid wasting power on outage events. These schemes are sub-optimal solutions of 



problem pi) . Based on the short-term schemes, two simple rules are discussed in the next 
subsections. 

1) Long-term truncated water-filling scheme: Similar to the short-term truncated water- 



filling scheme, we consider approximating Ix(p) in ( [37] ) by I tw (p) in ( f20| ), which results in 
the following problem 



Minimize 
Subject to 



(P(7)> 



p b > 0,6 = 1,..., B 



(58) 



Theorem 7: Problem ( [58] ) is solved by p(-f) given as 



Pb 



where rj is chosen such that 



mm 



lb' 

B 



1 

76 



B, 



Proof: See Appendix |VIII 



b=l 



(59) 



(60) 



Note that since / tw (p) is an upper bound on Ix(p), p(l) does not satisfy the rate constraint 
Ib(p(i),i) > R- By adjusting 77, we can obtain a sub-optimal p tw (-f) as follows 



Pb 



mm 



I 
lb' 



1 

lb 



1 B, 



(61) 



where now i] is chosen according to the true mutual information with discrete inputs, namely, 
we choose 77 such that 

1 B 

-Y^Ix(p t b w lb) = R- (62) 



6=1 
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Using this scheme, we obtain a power allocation p*™( 7 ), which is the long-term power 
allocation scheme corresponding to the sub-optimal p tw ( 7 ) of p opt ( 7 ). The performance 
of the scheme is illustrated by the dashed lines in Figures [5] and [6} As we observe, the 
performance of the truncated water-filling scheme is very close to that of the optimal scheme. 

Similar to the optimal power allocation scheme, we have the following duality between 
p£( 7 ) and p tw ( 7 ). 

Proposition 5: Consider transmission at rate R over the block-fading channel given in ([T]) 
with power allocation scheme pf™( 7 ) and long-term power constraint P(s). Then, indepen- 
dent of the fading statistics, the outage probability satisfies 

p out ( P \:,p(s),R) = p out ( P t s r( 7 ),s,i?), (63) 

where p*™( 7 ) is the truncated water- filling power allocation scheme with short-term constraint 

(p*r( 7 )> < s. 

Proof: See Appendix |IX] ■ 
Therefore, from Theorem [6] and Lemma |2} we have the following result. 
Corollary 2: Let pj^( 7 ) be the long-term power allocation scheme corresponding to £> tw ( 7 ) 
given in ( |6Tj ). Consider transmission at rate R over the block-fading channel given in ([T]) with 
the long-term power allocation scheme p{™( 7 ). Assume that the power fading gain 7 follows 
a Nakagami-m distribution given in ([5]). Then, the corresponding delay-limited capacity is 
non-zero if dp(R) > — , where 

R 



B 11- 



Proof: From Propositions [5] and |2} we have 



(64) 



PouMTh), P(s), R) = lC tw s-^ R \ (65) 

where d tw (R) > dp(R). Therefore, the proof follows as a result of Theorem [6] ■ 
In Figure [7] we also show in dashed lines the corresponding long-term truncated water- 
filling outage curves, and we observe the same asymptotic behavior as for the optimal scheme. 
2) Refinement of the long-term truncated water-filling scheme: In order to improve the 



performance of the sub-optimal scheme, we approximate Ix{p) by P (p) given in ( [30| ). 
Replacing Ix{p) in ( [37] ) by P ci (p), we have the following problem 

Minimize (p(l)) 

Subject to J Ef=i ^ ref ( Pblb) > R ( 66 ) 
p b >0,b=l,...,B. 
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Theorem 8: The problem given in <[66|) is solved by p(j) given as 



A 
ib ' 

K77, 

q 

lb ' 
77 - 

0. 



J_ 

lb ' 



n > 

I — K~f b 

-z- < v < — 

«±1 < „ < _9L 
lb — 1 — ^lb 

J_ < „ < schl 
ib — 1 ib 

otherwise, 



(67) 



where is chosen such that 



bib) 



BR. 



(68) 



Proof: See Appendix |Xj 
Following the arguments in Section 



6=1 



V-B.l 



we obtain the sub-optimal p mf (j) of p opt (l) 



from (67) by choosing r] in such that 



1 B 

6=1 



(69) 



The performance of the long-term power allocation corresponding to f/ ef (7), p[ t ef (7), is 
illustrated by the dashed-dotted lines in Figures [5] and [6] We observe that refined truncated 
water-filling leads to performance closer to that of the optimal schemes than truncated water- 
filling. The improvements are particularly clear for higher transmission rates. 

3) Approximation of Ix(p): The sub-optimal schemes in the previous sections are sig- 
nificantly less complex than the corresponding optimal schemes, while only suffering minor 
losses in outage performance. However, the sub-optimal schemes still require the computation 
or storage of Ix{p) f° r determining i]. This can be avoided by using an approximation of 
Ix(p)- Let Ix(p) be an approximation of Ix{p) and AR be the error measure given by 



AR = max\l x (p)-I x (p) 



(70) 



Then, for a sub-optimal scheme p('y), 77 is chosen such that 

B 

Y,ix{Pblb) = B(R + AR) 



b=i 



satisfies the rate constraint since 

B 



^{Pblb) > J2 ! ^(Pblb) - BAR = BR. 



(71) 



(72) 



6=1 
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Following [35], we propose the following approximation for Ix{p) 

I x (p) = M(l-e- c ^ 2 ) C3 . (73) 

For channels with QPSK input, using numerical optimization to minimize the mean-squared- 
error between Ix{p) and Ix{p), we obtain the parameters, ci, c 2 , c 3 , shown in Table [nj Using 
this approximation to evaluate rj in subsections V-B.l and V-B.2[ we arrive at computationally 
efficient power allocation schemes with little loss in outage performance. 

In Figure [8j we illustrate the significant gains achievable by the long-term schemes when 
compared to short-term schemes. As observed in [12], remarkable gains of 11 dB at an outage 
probability of 10~ 4 are possible with optimal long-term power allocation when compared to 
uniform power allocation for Gaussian input distributions. As shown in Figure [8j similar 
gains of the order of 12 dB at an outage probability of 1CT 4 are also achievable with discrete 
inputs. Note that, due to the Singleton bound, the slope of the QPSK-input short-term curves 
is not as steep as the slope of the corresponding Gaussian input or 16-QAM input curves. 
This is due to the fact that both Gaussian and 16-QAM inputs have SNR exponent d(R) = 4 
while QPSK has cLb(R) = 3. Figure [9] shows similar results comparing CM and BICM. In 
particular, the figure shows little loss between the corresponding power allocation schemes. 
This is due to the fact that the mutual information curves from CM and BICM with Gray 
mapping do not differ much [23]. Once again, in the case of BICM, the loss incurred by 
suboptimal schemes is negligible. 

We finally illustrate the application of the above results to practical OFDM channels. 
In particular, we show in Figure [10] the results corresponding to an OFDM channel with 
B = 64 sub-carriers, whose 9-tap symbol-period-sampled power delay profile is extracted 
from the ETSI BRAN-A model [36] using a zero-hold order filter. The power delay profile 



models a typical non-line-of-sight (NLOS) indoor office scenario and is given in Table III 



We observe a similar behavior as in the block-fading channel. In particular, we show that 
in practical OFDM scenarios impressive gains of more than 10 dB with respect to uniform 
power allocation (eventually reducing all outages) are possible. Note that due to the large 



frequency diversity induced by the time-domain channel in Table IIP the uniform and short- 
term power allocation curves do not reveal their respective asymptotic slopes in the error 
probability range shown in the figure. 
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VI. Conclusions 

We considered power allocation schemes for fixed-rate transmission over discrete-input 
block-fading channels with transmitter and receiver CSI under short- and long-term power 
constraints. We have studied optimal and low-complexity sub-optimal schemes. In particular, 
we have analyzed the optimal diversity orders and we have shown that, in the long-term 
case, outages can be removed provided that the short-term SNR exponent be greater than 
one. We have illustrated the corresponding performances, showing significant performance 
advantages on the order of 10 dB of the proposed long-term schemes when compared to 
uniform power allocation. Furthermore, we have shown that minimal loss is incurred when 
using the suggested sub-optimal schemes. We have also illustrated the applicability and 
performance advantages of the proposed techniques to practical OFDM situations. 
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Appendix I 

Optimal Power Allocation for Short-term Constraints 
Proof of Lemma Xrt Since p° pt (7) is the solution of ( [141 ), we nave 



B B 
6=1 b=l 

for any power allocation scheme p s t(7) satisfying the short-term power constraint. Therefore, 
Pr IxipTlb) < BRj < Pr ^ I x {p blb ) < BR^j , (75) 

and thus, 

Pont(pT(l), P, R) < Pout(Pb(l), P, R) (76) 

for any scheme p s t(7) satisfying the short-term power constraint. This proves that p° pt (7) 
is a solution of (fT3l). ■ 



Proof of Proposition [Tf With the optimal power allocation scheme p°t Pt (7)> the outage 
probability is given by 

Pont(pT(ll P, R) = Pr E < R^j ■ (7V) 

Since p°f(7) is the solution of ([14]), we have p° b pt > 0, b = 1, . . . , B and ± £]f =1 p° pt < P. 
Therefore, < p° b vt < BP, b = 1, . . . , B. Thus, Pout(p° pt (7)> Pi R) is lower bounded by 

Pout(p s T t (7) ) Pi R) > Pr f | J] ^(5?7i) < ^ (78) 

= Pou t (p cq (5P),P5,i?), (79) 

namely, the outage probability of block-fading channels corresponding to an equal allocation 
of power PB per block. Now, according to [27], under Nakagami-m fading statistics, we 
have that 

Pout(pT(7), Pi R) > K,(BP)- mdB{R) (80) 
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Conversely, since the power allocation scheme is optimal, the outage performance is upper 
bounded by the allocation scheme that assigns power P to each block. Therefore, 

PartfoTh), P, R) < Pout(p cq (P), P, R) (82) 

= JCp- md B(R). (83) 



From ( |8T| ) and ( [83] ), we have 

Pout(pTh),PR) = iCo P tP- mdB{R) . 



(84) 



Thus, the diversity obtained by the optimal power allocation scheme is given by mc?s(i2), 
which is the same as that of the uniform power allocation scheme. This concludes the proof 
of the Proposition. ■ 
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Appendix II 
Proof of Theorem[2] 

Proof: The power allocation algorithm of interest is the solution of the optimization 



problem d2Tl). Since f(pblb) is constant at log 2 (l + (3) for p b > having pb > — does 

■ ' lb lb 



not provides any gain to the target function in pi) . Therefore, the solution of the following 
optimization problem 

Minimize /o (p) - - i E&Li 1o S2 ( 1 + Pblb) 
Subject to fb{p) = ~Pb < 0, 



9b{P) 



Pb 



76 



<o, 



1,...,B 
1 B 



(85) 



Hp) = EtiPb<BP 



is also a solution of (21 ) 



It can be verified that the functions /o(p)? h(p), 9b(p), b = 1, . . . , B, and h(p) are convex. 
Therefore, according to the Karush-Kuhn- Tucker (KKT) conditions [37], the solution p tw of 



(85) must satisfies 



v > 
BP I = 



.6=1 



=i 

B 

^2pT-bp<o 



b=i 



lb log 2 e 



pT > o, 


6 


= 1,- 


..,B 


x b > o, 




= 1,- 


..,B 


-^pT = o, 




= 1,. 


..,B 


a b > 0, 




= 1,. 


..,B 


pT - - < o, 




= 1,. 


..,B 


( rfw -f)^ - 




= 1,. 


..,B 


A b + a b + z/ = 0, 


b~- 




.,B, 



(86) 
(87) 

(88) 

(89) 
(90) 
(91) 
(92) 

(93) 
(94) 
(95) 



l+pflb 

where z/, a b ,b = 1, . . . , B are the Lagrangian multipliers. For any b, 

. If X b > 0, from ([91), pf = 0. Therefore, from ([94]) a b = 0. In this case condition ([95]) 



is satisfied only if v > % log 2 e. 
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• If Xb = 0, we have the following cases 

- If a b > 0, from ( |94| >, p£ w = K In this case, (|95]) is satisfied only when v < lbl °^l e . 

- If a b = 0, from <|93), p£ w = ^ - ±. Furthermore, from (f89j) and ([93]), we have 

<^<76log 2 e. 
Therefore, for any choice of v, we must have 



v < 



log2 e 1 76 ' 



mm 



76 ' /3+1 

otherwise 

/ log 2 e _ _1_ 
V V ^ 76/ + 



76 log 2 e 
/3+1 

— < f < 7fo log2 e 



(96) 



(97) 



The solution in (]97[) satisfies conditions (|89|) — (|95|). We are left to choose v > such that 

, P. Therefore, 



conditions ( |86| ) - ( [88] ) are satisfied. If v = 0, from ( |97| ), p£ 
z/ = is valid only if 



tw _ £ b 

76 ' 



1. 



B 



(98) 



Otherwise, if 5^&Li ~ > -^-P> choose ^ such that 



ib 



B 



£ 

6=1 



= PP. 



(99) 



Therefore, letting rj 

Ef i A < PP, 



Otherwise, the solution is 



where rj is the solution of 



logo e 



the solution of (85) can be summarized as follows. If 



7fe 



P. 



tw 

p b = mm 



76 ' 



lb 



(100) 



(101) 



E 

6=1 



mm 



lb 



V - 



lb 



BP. 



(102) 
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Appendix III 
Proof of Proposition [2] 

Proof: With the truncated water- filling power allocation scheme p*™(7), the outage 
probability is given by 



P out (p t s r(7), P, R) = Pr ^ J2 J ^(P6 W 7 6 ) < rJ 



(103) 



The outage probability can be upper bounded by 



P out (p^(7), P, R) < Pr (± ^PTlb) < Rj , 



(104) 



where I^{p) is a lower bound to given by 



4(p) 



0, otherwise. 



(105) 



We now further upper bound P out (p*™(7), P, R) using the following proposition. 



Proposition 6: Consider the truncated water- filling scheme given in ( |22[ ). For any channel 
realization 7, we have 

B B 



6=1 



6=1 



where I^(p) is given in ( |105| ). 

Proof: According to ( |105[ ), I^(P / -fb) is non-zero only if 7^ > -|. Therefore, we need 
to prove that if 7^ > p then > for all realization of 7. Without loss of generality, 
assume that 71 < . . . < 75. If 7b < 75, ( |106[ ) is certainly true. Otherwise, there exists a 
k,l < k < B, such that 7^.1 < 4 < 7fc < . . . < 7s- Consider the following two cases: 



• If Ef=i i < x then from ©' ^ w = - 



76 ' 



1 P. 



Otherwise, from d22|), the power allocation solution is given by 



mm 



I 
lb 



V - 



1 

lb 



B, 



where r] is chosen such that 



Erf 

6=1 



tw 



PP. 



(107) 



(108) 



Since 7^ > 75, we have from ( |107[ ) 



pf<^<^<P, b = k,...,B. 
7fe 7fe 



(109) 
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Therefore from (11081), 



6=1 



Now, suppose 



then, for b — 1, . . . , k, 



b=k+l 

(3 + 1 



T] < 



Ik 



tw . 1/3 + 1 1 
P b <V — < 







< P. 



(110) 



(111) 



(112) 



lb Ik Ik Ik 
Thus, J2b=i Pi™ < kP, which contradicts to ( | 1 10| ). Therefore, assumption < | 1 1 1 1 ) is invalid. 
We then conclude that r) > ^ > @±±,b = k, . . . , B. Therefore from ( [107] ), p£ w 

£ 6 = k B. 



ib ■ 



Therefore, in all cases, we have p^ = ^ if ^ b > J^. This concludes the proof of the 



proposition. 

From Proposition [6j we can further upper bound P out (pgf (7), P, R) by 



P out (p t s r( 7 ), P, R) < PLMP), P, R) = Pr (± J2 T x( p ^ < R j ■ 



(113) 



The asymptotic behavior of P^ ut (p CC[ (P) , P, R) is given by the following proposition 
Proposition 7: Assume that 7 b follows the distribution given in (|5]), then P ( f ut (p eq (P), P, R) 
in ( |1 13[ ) asymptotically behaves as 

Pl( Pcq (P),P,P) = /C /3 p- m ^, 



where 



d g = 1 + 



P I 



P 



(114) 



(115) 



and Ix (p) is the input-output mutual information of a AWGN channel with input constellation 
X and SNR p. 

Proof: Consider the random set given by S s = {« e {1,...,P} : 7^ > -|}. Then for 



6=1, 



P, 



Pr(6 e S 3 ) = Pr 7fc > 



P 



1 — Pv 



(116) 
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The asymptotic behavior of pp is given by 



T ( m, m— 



r(m) 

r(m) - i (mf )* 



1 -P/S 



T(m) 



p-m 



Since 71, ... ,7b are independent random variables, jS^I is binomially distributed 

'P N 



Pr(|^|=t) 



,B-i 



t ) \ T(m) 



p-m(B-t) 



Now from (105), 



0, otherwise. 



Therefore, 



P out (P,R) = Pr f l5(P 76 ) < BR 

= Pt(\S p \I x (J3)<BR) 

, , . „ ■ BR 

Pr 5 fl < 



£ Pr (l^l=t) 



r br 
I 

£ 



-1 



p-m(B-i) 



At high P, the dominating term in ( |127| ) is the term with t = t± 

P out (P,R) = ICpP- md ^ R \ 

where 



BR 



B 1 



P 



^(P) = P - ti = 1 + 

This concludes the proof of the proposition. 
Finally, from ( |1 13| ) and Proposition [7J we have 

P out (p^(7), J P,^)</C /3 P- m ^^ 

as required by the Proposition. 



W)J\ ' 



(117) 

(118) 
(119) 

(120) 
(121) 

(122) 

(123) 

(124) 
(125) 

(126) 

(127) 

1. Therefore, 
(128) 

(129) 



(130) 
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Appendix IV 
Proof of Theorem[3] 

Proof: Similar to Theorem [2| a solution Pst f (7) to the optimization problem given in 
3Tb satisfies the KKT conditions [37] for the following problem: 



Minimize - Y. b =i ^ (Pblb) 
Subject to J2b=iPb — BP 

p b <f b ,b = l,... 

Pb > 0,b = 1,. . . . 

v > 

B 

J2pf < BP 

6=1 

- BP 



Therefore, p ref (7) satisfies 




X b > 0, 

pf > o, 



rcf 



n l P f 



X bP v b 
n > o, 

Pf < K 
lb 

lb 



and 



76 log 2 g 

Pl Ci 



- \ b + r b + v 
Xb + n + v 



o. 
o. 



1J ^> -Xb + n + u >^^. 

l+Q tJ a 1 

for 6 = 1,..., £. 

For any 6, consider the following cases 
. If X b > 0, then from ([137]), ( [UP] ), we have pf 

is satisfied only if v > j b log 2 e. 
. If A 6 
. If X b 



B 



1,...,B 
1,...,B 
1,...,B 
1,...,B 
1,...,B 

1 B 



if pf < 

<*. < 
it 

if pf 



76 

if a < ref < A 

76 ^" — 76 



7b 



(131) 



(132) 
(133) 

(134) 

(135) 
(136) 
(137) 
(138) 
(139) 

(140) 



(141) 



r b = 0. In this case, condition (141 ) 



and r b > 0, then from ( |140| ), p T b ' 
n 



rcf 



^76 log 2 e 



E -. Therefore, from ( 141 ), v < 

76 ' " < P 



0, from ( 141 ), we have 
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+ Pf 



v 76 
?2 r 



- when < pf < — or equivalently when 7 '' 1 1 °^ 2 6 < v < 75 log 2 e 



+ pf = if ^ < K ef < £ ^ 



7b 

g log2 e 7ft < < «76 log2 e 



+ Pb 



ref _ y 7ft K l°g2 e 



< 1/ < 



7ft 'Og2 e 

7d a — " — 1+a 

Therefore, for any choice of v, we have 

A 

7b ' 



Pf 



V < 



«7i> l°g2 e 



/3 



* log2 e 



7ft ' 



«7ft l°g2 e 
/3 

«7ft log2 e 



< 1/ < 



< V < 



K~jb log2 e 
a 

7ft log2 e 



log2 e L 7ft log2 e 

1/ 7 b ' 



0. 



1+a 

1+a < f < lb log 2 e 
otherwise. 



(142) 



We are left to choose v > such that conditions ( |132[ )-( (T34| ) are satisfied. If v = 0, then 
from ( |142[ ), pf = — , 6 = 1, . . . , B. Furthermore, from ( |133[ ), v = is valid only if 



B 







^ ib 



(143) 



6=1 



If Ef=i - > Bp , t nen v > 0. Therefore, from ( fT34l >, z/ is chosen such that 

B 

J2pf = BP. 

6=1 

Therefore, by denoting r\ = we obtain p rcf (7) as defined in the Theorem. 



(144) 
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Appendix V 
Proof of Theorem [4] 

We first consider the following Proposition, which is a generalization of the result in [12] 
to channels with discrete inputs. 

Proposition 8: The solution of ( [35] ) has the following form 



PT(7) 



p opt (~f), with probability w("y) 



(145) 



0, with probability 1 — $(7), 

where p opt is the solution to the problem in ( [37] ) and w(j) is the solution of 

Maximize E [10(7)] 

Subject to < ^(7) < 1 

E[(p^(j))w(j)]<P. 
Proof: From ( |145[ ) and ( |146[ ), we have 



(146) 



E [<Pit Pt (7)>] = E [(p opt (7))^(7)] <P, 



(147) 



which shows that p° pt (7) satisfies the long-term power constraint. We need to prove that 

^out(pr t pt (7), P, R) = 1 - E H 7 )] < P out (p( 7 ), P, R), (148) 

where p(7) is an arbitrary power allocation scheme satisfying the long-term power constraint 

E[(p( 7 ))]<P. 

Given a channel realization 7, define the region 



1 B 1 

£5>(W6)>i? 
6=1 J 



(149) 



and 



W ( 7 )4p r (p( 7 ) G*4( 7 ,P)). (150) 

Since .4.(7, P) is the power allocation region that does not cause outages, the outage proba- 
bility given a channel realization 7 is 1 — u>(/y), and the overall outage probability is given 
by 

P out (p(7),P,i?) = l-EH7)]. (151) 



We now prove that w(y) satisfies the constraints of the problem given in ( |146[ ). By definition 



( |150| ) we have that < w(j) < 1. Furthermore, since p opt (7) is a solution to ( [37] ), we have 

V p( 7 ) G Afr, R), (fP opt (7)> < (P(7)>- d52) 
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Therefore, conditioned on 7 , the expectation of (p(7)) over the distribution of p(7) can be 
lower bounded as follows. 

Ep (7)eR B [<p(7)}|7] > Ep(7)^(7.«) [<P(7))l7] (153) 
>(fp opt ( 7 ))Pr(p( 7 )e.4( 7 ,P)) (154) 
= (^ opt ( 7 )M 7 ). (155) 



Thus, since E [(p(7))] = E 7eK? E p(7)eM s [(p( 7 ))| 7 ] 



< P, we have 



E [(^(7)^(7)] < P. (156) 



As a result, w("y) satisfies the constraints in ( |146| ), and thus, 

E [w{j)} > E [w{j)] . 

Therefore, we finally have 

P out (p° pt , P, R) = 1 - E [w(j)] < 1 - E [w( 7 )] = P out (p, P, P) (157) 



for any arbitrary p(7), which shows that p° t pt ( 7 ) is a solution to the problem given in ( [35] ) 



We have the following Proposition, which gives the solution to the problem given in ( 146) 



Proposition 9: Suppose that % follows a continuous probability density function. Then, a 
solution to the problem in ( |146[ ) is given by 



1, if (p opt (j)) < s 
"■(7)=< d58) 
0, otherwise, 



where s satisfies 

I s = 00. if lim^^ P(s) < P 

(159) 



P(s) = P, otherwise, 
and 

P( S )4E[(^( 7 ))%)]. 
Proof: If lim^oo P(s) < P, $(7) = 1 (which corresponds to s = 00) is certainly a 

solution to the problem. 

Consider the case when lim^oo P(s) > P. We first prove the existence of an s satisfying 



( |159| ). Denoting ^7(7) as the pdf of 7 , we can write P(s) as 

P(s) = [ (p opt (7))/ 7 (7)^7, (160) 
Jn(s) 
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where TZ(s) is defined in p9y . For all s > s,1Z(s) C 1Z(s ). Therefore, from ( |160[ ), P(s) is 
an increasing function of s. Due to the continuity of the fading statistics and of the mutual 
information curve, (p opt (j)) is a continuous function of 7. 

Without loss of generality, assume that 71 > . . . > j B . We first prove that p° pt > 0. This 
is in fact the case since for any power allocation scheme p(j) such that pi (7) = 0, and 
Pk{l) > satisfying the rate constraint, the power allocation scheme p'ij) with 



P*(7)£: 



b = 1 

Phil) = \ 0, b = k (161) 

Pb(l), otherwise, 

which has < also satisfies the rate constraints. 

Now assume that 7' satisfies 7^ > 71, 7^ = 7&, b = 2, . . . , B. Consider the power allocation 
scheme p(j') satisfying 



Pb{l') 



7i 



(162) 



P°b (7), 



otherwise. 



Obviously, 



6=1 6=1 

Therefore, (£> opt (7')> < (p(l')) < (p opt (l)) (since p° pt (7) > 0). This proves that (p opt (j)) 
is a strictly decreasing function of 71 for any fixed 72, . . . , 75. 

Due to the aforementioned monotonity and continuity of (p opt (j)), given 72, • • • , 7s, there 
exists a unique 71 (s, 72, • • • , 7b) such that (p opt (7)) = s for any s > 0, and we can rewrite 
the region TZ(s) as 

1l{s) = {7e»f -Ix > 7i (s, 72,..., 7b)}- (163) 



Additionally, for all e, there exists a 5(e) such that 

71 (s + 6,72, ... ,7s) = 7i0,72, 
linage) = 



. 1b) + 6{e) 



eS(e) < 0. 



(164) 
(165) 
(166) 
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Therefore, denoting 7 S = 71 (s, 72, • • • , 773), we have 

limP(s + e) = lim / (p opt (l)) f 7 h)dj (167) 

= lim / ( ^ (p opt h))hh)d7i) d l2 ... d lB (168) 

= / f lim T (^ opt ( 7 ))/ 7 (7)^7iV^---^ ( 169 ) 
= I ( r(p° pt h))hh)dli) d l2 ... d lB (170) 

^72, ■■■,7s V</7s / 

= P(s). (171) 

Thus, -P(s) is an continuously increasing function of s, which proves that there exists an s 
satisfying P(s) = P since lim^oo P(s) > P. 
On the other hand, for any 10(7), we have 

E [(^(7)^(7)] - P = E [(p opt ( 7 )M7)] - E [(^(7)^(7)] (172) 

= / W ( 7 )(p opt (7))^ 7 (7) (173) 

^Mf\7e(s) 

-/ (l- W ( 7 ))(p opt (7))^ 7 (7) (174) 
Jn{s) 

w( 1 )dF 1 ( 1 ) - [ (1 - w( 1 ))dF 1 ( 7 ) ) (175) 

Ms) J 



> s 



\Tl(s) 



s(E[w( 7 )]-E[w( 7 )}), (176) 



where (173) and (176) are due to 



1, if 7 e TZ(s) 
u>( 7 ) = { (177) 

0, otherwise, 



and ( |175[ ) is obtained using the following bounds 

(p opt h)) <s, if 7 G (178) 

(£, opt (7)> > s, if 7 ^ (179) 

Therefore, for all 10(7), E [117(7)] > E [^(7)] implies E [(p opt (7)M7)] > P, which 
violates the problem constraint. Thus, u?( 7 ) is a solution to the problem. This concludes the 

proof of the proposition. ■ 

The proof of the Theorem is obtained by applying Propositions [8] and [9j 
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Appendix VI 
Proof of Theorem [5] 

Proof: Since (p), — J2b=i ^x(plb), — Pb are convex functions of p(j), applying the 
KKT conditions [37] and note the fact that [22] 



4-Ixip) = -^MMSE*(p), 
dp log 2 

the solution p opt (i) of ((37]) satisfies the following conditions 

v > 

B 



B 



(180) 

(181) 
(182) 

(183) 



6=1 



\ b > o, 

opt < 0j 
opt __ Q 



^ 7b MMSE*(p°%)-A fe = 



log 2 



where u,X b ,b = 1, . . . , B are the Lagrangian multipliers. Letting r\ 
have 



(184) 

I B (185) 

(186) 
B (187) 
, for any b, we 



1) • • • ; 

v 

log2 



• If Xb > 0, then from ( |186[ ), p^ pt = and thus ( |187| ) requires 77 < -i. 



. If A 6 = 0, then from (187), 



opt 



lb 



-MMSE7, 1 



1 



(188) 



and rj > ^ since MMSE^p) < 1. 



Therefore, with any choice of 77, we have 



, °Pt 

Pb 



-MMSE^ 1 , 7fe > i 



lb 




76 



-MMSEJ, 1 [ min < 1, 



otherwise 
1 

Vlb 



(189) 



(190) 



for 6 = 1, . . . , B. We are left to choose 77 > such that (182) and (183) are satisfied. From 



(189), (183) is not satisfied if 77 = 0. Therefore, from (182), 77 is chosen such that 



B 



(191) 



6=1 



as required by the Theorem. 



February 5, 2008 



DRAFT 



Appendix VII 

Asymptotic Analysis of Power Allocation for Long-Term Constraints 

Proof of Proposition [3} From the definition of differentiation, we have 

d , N , P(as) - Pis) 

—P{s) = lim ^, (192) 

as all as — s 

where 

P(as) = E Jen{as) [<p( 7 )>] (193) 
= E 7e ^ {s) [<p( 7 )>] + E 767 * MW [(p( 7 )>] (194) 
= P(s) + E imasm(s) [(p( 7 )>] . (195) 

Note that V 7 G TZ(s), we have that (p( 7 )) < s. Therefore, since a > 1, (p( 7 )) < as, which 
implies that 7 G TZ(as) and thus, 72(s) C TZ(as). Now, let / 7 ( 7 ) be the pdf of the 7 . Since 

V 7 G lZ(as), (p>( 7 )) < as, we have 

E 7 ^(a S )\^(s) [<p( 7 )>] = / {P{l))h{l)dl (196) 

^ 7 e7e(as)\7e(s) 



< as / / 7 ( 7 )d 7 (197) 

</ 7 e7e(as)\7e(s) 

= as [Pr ( 7 i K{s)) - Pr ( 7 £ 72(as))] . (198) 



From the assumption in ( |48| ), and noting that P out (p lt ( 7 ), P(s), 7?) = Pr( 7 ^ 72(s)) = 
)Cs~ d ( R \ we have 

E 7e7 * (as ) Ws ) [<P(7)>] < osK {s~ d{R) - (as)- d ^) . (199) 
On the other hand, since V 7 G TZ(as) \ 72(s), (p(j)) > s, by similar arguments, 

E le n(as)\ns) [<p(7)>] > s£ - (as)- d ^) . (200) 



Therefore, from ( fl92l , ( [195) , ( fT99j ), ( |2Q0] >, we have 



lim ^ ^ < — P(s) < lim ^ '-. (201) 

0J.1 a — 1 as oil a — 1 



Since 



lim ^ ^ = lim ^ '- = Kd(R)s- m , (202) 

ail a — 1 ail a — 1 

we have that 

■f-P(s) = JCd(R)s- d{R) (203) 
as 

which concludes the proof. ■ 
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Proof of Theorem pr We begin with the first part of the theorem, i.e. lim^oo P(s) = 
Pth < oo when d(R) > 1. From Proposition [3j we have 

lim ( ^P(s)] s d{R) = fCd(R). (204) 

s^oo yds J 

Therefore, for any e > 0, there exists a finite s\ such that for all s > s%, 

■f-P(s) ] s d{R) < Kd(R) + e, (205) 
ds I 



or equivalently, for all s > Si, 



Thus, 



-f P(s) < (Kd(R) + e)s- d(fi) . (206) 



lim P(s) = P(si) + lim / f ^-P(t) ) dt (207) 
s^oo s^oo y ydt J 

< P( Sl ) + lim / (K,d(R) + e)t- d{R) dt (208) 
s ^°° J Sl 



P(*i + lim ^ — — '-, (209) 

s^oo l — d{R) 



which gives 



= Pth<oo, (211) 
when d(P) > 1 as required. Furthermore, from the definition of the long-term exponent, 

r ,px -\QgP out (Plt(l),P(s),R) 
MR) = P(^oo \^P{s) (212) 

lim -togP^^PW^) (213) 
log P(s) 



Therefore 



if d(P) > 1. 



-log (ICs-^) 
lim 6 1 ^ (214) 

s^oo log P(SJ 

r d(P)logs 

hm — — . (215) 

s^oo logP(s) 



d lt (R) = hm rf(i?) D lQg(g) = oo (216) 
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In the second part of the theorem, where we have d(R) < 1, then from ( |209| ) we observe 
that linis^oo P(s) = oo. Applying L'Hopital's rule to ( |215[ ), we obtain 

d lt (R) = lim (217) 

Applying Proposition [3j we can further write d\ t (R) as 

d "W = ,^OT)- (218) 

Further applying L'Hopital's rule and Propostion [3] yields 

^ = a c(i-S'U (219) 



(220) 



1 - ' 

which completes the proof. ■ 
Proof of Proposition^ From [21], the optimal power allocation scheme for system with 
short-term power constraint s is given by 



opt 

Pb 



(7) = ^MMSE^ 1 (min jl, J) , b = 1, . . . , B, (221) 



where r\ is chosen such that the power constraint is satisfied, 



B 



Y,PT(l) = Bs. (222) 

6=1 

Transmission with rate R and power allocation scheme p°t Pt (7) is in outage if and only if 
there is no rj satisfying 

Ix (MMSE^, 1 ^min j 1, ^ > BR (223) 

^ -^MMSE^ 1 ^min |l, < Bs. ( 224 ) 

Similarly, from Lemma [5] and Theorem |4j transmission with power allocation scheme p° t pt (7) 
is also in outage if and only if there is no r\ satisfying ( |223[ ) and ( |224[ ) simultaneously. 
Therefore, 

P out ( Plt ( 7 ), P(s), R) = Pout( P r (7), s, R) (225) 
as required by the Proposition. ■ 
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Appendix VIII 
Proof of Theorem [7] 



Proof: Similar to Appendix |n} a solution to the problem given in ( [58] ) is given by 
solving 

Minimize (p(-y)) 

Subject to J2b=i lo g2(! + Pblb) > BR 



Pb 



< 



(226) 



p b >0, b=l,...,B. 



The problem given in (|226J is a standard convex optimization problem. Therefore, according 
to the KKT conditions [37], a solution p*(~f) to the problem satisfies 



v > 

BR-J2^og 2 (l + pH)<0 



b=i 



v 



BR- ^log 2 (l + p£ 76 ) j =0 

v 6=1 / 



(227) 
(228) 

(229) 



A 6 > 0, 


6 = 1,. 


.,B 


(230) 


P*b>0, 


6 = 1,. 


.,B 


(231) 


hpt = °) 


6 = 1,. 


.,B 


(232) 


«6 > o, 


6 = 1,. 


.,B 


(233) 


pJ--<o, 

7ft 


6 = 1,. 


.,B 


(234) 


a(p*-^) = 0, 
7b 


6 = 1,. 


.,B 


(235) 


7fclog 2 e 
i/ A fo + a fe = 0, 


6 = 1,. 


.,B 


(236) 



where u, A&, otb, b = 1, . . . , B are the Lagrangian multipliers. For any 6, 

0. Therefore, (|236[) requires 



• If At > 0, from ( |232[ ) and ( |235| ), we have p£ 



Ob 



V < 



1 



76 l°g2 e ' 



• If Afe = 0, a b > 0, from ( |235[ ), we have p£ 



76 



Therefore, (236) requires z/ > 



76 log 2 e ' 



• If A;, = «6 = 0, from ( |236[ ), we have pi = u\og 2 e — Therefore, ( |230[ ) and ( |234[ ) 



require — A — < v < ^ +1 . 

n 76 l°g 2 e — — 76 l°g2 e 



February 5, 2008 



DRAFT 



43 



Therefore, for any choice of v, we have 



Pb 



76 ' 

v log 2 e 
0. 



v > 



/HI 



76 log 2 e 
1 < V < 



mm 



I 
lb' 



j_ 

76 ' 76 log 2 e 

otherwise 
1 



76 log 2 e 



v log 2 e 



76 



(237) 



(238) 



We are left to choose v such that the conditions in (227)-(229) are satisfied. From (238), 



( |228| ) is not satisfied if v = 0. Therefore, from (J229J), we choose v such that 

B 



B 



lb 



v log 2 e - 



76 



BR. 



6=1 6=1 

Finally, denoting 7] = z/log 2 e, we have p*(7) as required by the Theorem. 



(239) 
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Appendix IX 
Proof of Proposition [5] 

Proof: According to Theorem [2j the truncated water-filling scheme with short-term 
power constraint s can be written as follows 



P 



/fast, 7b) 



with 77^ chosen such that 




where 



/fa, 7) 



oo if y;f 1 - < Bs 

= Bs otherwise, 



( 1 

mm < — , 7] 

.76 V 7fe 

B 



(240) 

(241) 

(242) 
(243) 



Similarly, from ( |6T| ), <p% 



g(v,nr) = $^/fa,7&)- 

6=1 

/fait, 7b), & = 1, . . . , B with % chosen such that 

7fa lt , 7 ) = BR, (244) 

where /fa, 7) = X)&Li £tr(/fa, 7&)7&)- Transmission with power allocation scheme pf t w (7) is 
in outage if #fait,7&) > Bs. 

Consider truncated water-filling schemes with (3 chosen such that Ix(P) > R. Consider a 
channel realization 7, and, without loss of generality, suppose 71 > . . . > jb. If transmission 
with power allocation scheme p*f(7) is in outage then 



Vst < 



(3+1 



IB 

9 fast, 7) = Bs 
-ffast, 7) < BR. 

Noting that /fa, 7) and gfa,7) are increasing function of r] for rj < ^±r, from (|244j) and 



(245) 

(246) 
(247) 



JB 



< |247| ), we have % > r] st and thus, gfait,7) > gfa s t,7) = Bs. Therefore, transmission with 
power allocation scheme p\™(j) is also in outage. 

By similar arguments, we also conclude that if transmission with p{^(7) results in an 
outage event, then transmission with p*™(7) also results in outage. 

Therefore, 

P out (pf t w ( 7 ), P(s), R) = P ou t(p*r(7), s, R). (248) 
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Appendix X 
Proof of Theorem [8] 

Proof: Similar to Theorem [2| a solution p*(*y) to problem ( |66] ) satisfies the KKT 
conditions [37] for the following problem 

Minimize 

Subject to Ef=J mf (Pblb) > BR 
Pb<£, b = l,...,B 
p b >0, 6 = !,...,£. 



Therefore, satisfies 

v > 

B 

BR-J2r cf (pH)<o 

6=1 



6=1 



X b > 0, 
Pt>0, 

hpt = o, 

pt--<o, 
lb 



n\pb- 



ft_ 
lb 



o. 



and 



1 - v 



lb log2 e 
, K log 2 e 



- Afe + r b = 0, 



l-z/™-A b + r 6 = 0, 



7b log2 e 
1+T 



> 1 - X b + r fe > z/ 



7 b Klog 2 e 



for 6 = 1,...,£. 

For any 6, consider the following cases: 



6=1,.. 

if < f b 
if t < ^ < 
if pI = 1 



B 
B 
B 
B 
B 

B 



(249) 



(250) 
(251) 

(252) 

(253) 
(254) 
(255) 
(256) 
(257) 

(258) 



(259) 



If Xb > 0, from ( |255| ) and ( |258| ), we have pi = n = 0. In this case, condition (|259j is 



7t l°g 2 e 



satisfied only if z/ < 
If A 6 = and r b > 0, from ( gggfr , pt = - Therefore, from ([259]), i/ > . 
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• If Xb = r b = 0, from (12591), we have 



+ pi = v log 2 e — - when < p\ < — or equivalently, when 



+ pi = z/Klog 2 e when ^ < pi < ^, or equivalently when K7fcl ° g e 
+ p £ = £L when A tB ~. < v < 



— when , 

7b 7b i°g 2 e 



< i/ < 
< v < 



a+l 



Klb log 2 e " 



«7b !og2 e ' 

Therefore, for any choice of z/, we have 



Pb 



i/k log 2 e, 

< -, 
v log 2 e - 
0. 



v > 



«7b log 2 e 



«7b log 2 e 

1+q 
7b log 2 e 



< 1/ < 



< 1/ < 



1 



1 



7b ' 7b l°g 2 e 

otherwise. 



< i/ < 



K 7b l°g 2 ' 



«7b log 2 e 

q+1 
7b log 2 e 



(260) 



We are left to choose v such that conditions ( |250| )— ( |252| ) are satisfied. From (J259J), p£ = 
0, b = 1, . . . , B if f = 0. Therefore, ( |251[ ) requires that v > 0. Thus, from ( |252[ ), we need 
to choose i/ such that 

Y,I rc \ptlb) = BR. (261) 



6=1 



Therefore, by denoting r\ = ulog 2 e, we obtained as defined in the Theorem. 
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TABLE I 

Parameters po,n,a and a for the refined power allocation scheme. 



Modulation Scheme 





QPSK 


8- 


PSK 


16-QAM 


64-QAM 




CM 


BICM 


CM 


BICM 


CM 


BICM 


CM 


BICM 


Po 


3 


3 


7 


7 


15 


15 


63 


63 


K 


0.3528 


0.3528 


0.4693 


0.4744 


0.56 


0.5608 


0.6581 


0.6460 


a 


1.1327 


1.1327 


1.1397 


1.1234 


1.347 


1.3452 


1.5255 


1.5978 


a 


1.585 


1.585 


2.1677 


2.0922 


5.8884 


5.8264 


18.954 


19.8884 
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TABLE II 

Optimized a, c 2 and c 3 parameters for the approximation (73) of [35]. 



Modulation Scheme 



QPSK 8-PSK 16-QAM 64-QAM 





CM 


BICM 


CM 


BICM 


CM 


BICM 


CM 


BICM 


Cl 


0.77 


0.77 


0.61 


0.81 


0.48 


0.59 


0.47 


0.4 


C2 


0.87 


0.87 


0.68 


0.06 


0.61 


0.06 


0.44 


0.05 


C3 


1.16 


1.16 


1.45 


1.75 


1.48 


1.65 


1.87 


1.63 


Ai? 


0.0033 


0.0033 


0.0241 


0.0223 


0.0414 


0.0259 


0.0977 


0.0656 
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TABLE III 

Power delay profile of the normalized ETSI BRAN-A channel model using a zero-hold order filter. 



Delay (multiples of 50ns) 


Normalized path power (dB) 


1 


-3.4630 


2 




3 


-8.9151 


4 


-12.8223 


5 


-19.9222 


6 


-21.1202 


7 


-25.4329 


8 


-29.7891 


9 


-34.1993 
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Fig. 1. Block diagram corresponding to the channel system model with CSI at the transmitter and the receiver. 
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Fig. 2. Mutual information of QPSK, and the approximations used by truncated water-filling and its corresponding 
refinement. 
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P(dB) 

Fig. 3. Outage performance of various short-term power allocation schemes for QPSK-input block-fading channels with 
B = 4 and Rayleigh fading. The solid-lines represent the optimal scheme; the solid lines with o represent uniform power 
allocation; the dashed lines and dashed-dotted lines represent truncated water-filling and its corresponding refinement, 
respectively; the dotted lines represent the classical water-filling scheme. 
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5 10 15 20 



P (dB) 

Fig. 4. Outage performance of various short-term power allocation schemes for QPSK-input block-fading channels with 
B — 4 and Rayleigh fading. The solid-lines represent the optimal scheme; the solid lines with o represent the uniform power 
allocation; the dashed lines and dashed-dotted lines correspondingly represent the truncated water-filling and its refinement 
with j3 = 15. 
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-6 -4 -2 2 4 6 8 10 12 14 

P(dB) 



Fig. 5. Outage performance of various long-term power allocation schemes for QPSK-input block-fading channels with 
B — 4 and Rayleigh fading. The solid-lines represent the optimal scheme; the dashed lines and dashed-dotted lines 
correspondingly represent the long-term truncated water-filling (pft"(7) with Pi) and its refinement (p[ t ef (7) with fi-i). 
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Fig. 6. Outage performance of various long-term power allocation schemes for 16-QAM-input 4-block block-fading 
channels under Rayleigh fading. The solid-lines represent the optimal scheme; the dashed lines and dashed-dotted lines 
correspondingly represent the long-term truncated water-filling (pft"(7) with Pi) and its refinement (p[ t ef (7) with j3i). 
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P(s),s(dB) 

Fig. 7. Outage performance of long-term power allocation schemes as a function of P(s) and s (dual short term scheme) 
with QPSK-inputs in a block-fading channel with B = 4, m = 0.5, R = 1.7 and Rayleigh fading where cIb(R) < ^. 
Solid-lines correspond to the optimal scheme and dashed lines correspond to truncated water-filling with /?i = 3. The curves 
as a function of P(s) are the leftmost. 
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P(dB) 



Fig. 8. Outage performance of various short- and long-term power allocation schemes in a block-fading channel with 
B — 4, R — 1, Rayleigh fading and Gaussian, QPSK and 16-QAM inputs. The thick solid line corresponds to the Gaussian 
input; the thin solid-lines represent optimal scheme; the solid lines with o represent uniform power allocation; the dashed 
lines and dashed-dotted lines represent truncated water-filling (pfi™(7) with f3\ = 6 dB for QPSK input and 14.3 dB for 
QAM input respectively) and the dashed dotted lines represents the corresponding refinement (pf t of (7) for QPSK input with 
#2 = 5.5). 
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16-QAM 



16-QAM 
Gray. 



16-QAM. 



16-QAM 
V Gray 



Long-term 
L 



Short-term 



-5 



5 

P (dB) 



10 



15 



Fig. 9. Outage performance of various short- and long-term power allocation schemes in a block-fading channel with B = 4, 
R = 1, Rayleigh fading, Gaussian and 16-QAM CM and BICM (with Gray mapping). The thick solid lines correspond to 
the Gaussian input; the thin solid-lines represent optimal scheme; the dashed lines represent truncated water-filling (pft"(7) 
with /3 = 13 dB). 
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Long-term 



Short-term* 



-5 



5 

P(dB) 



10 



15 



Fig. 10. Outage performance of various short- and long-term power allocation schemes in an OFDM channels with 
B — 64 carriers, R = 1, Rayleigh fading and Gaussian, QPSK inputs. The thick solid line corresponds to the Gaussian 
input, the thin solid lines represent the optimal scheme; the solid lines with o represent uniform power allocation; the dotted 
line represents the pure water-filling; the dashed lines and dashed-dotted lines respectively represent long-term truncated 
water-filling pf^(7)and its corresponding refinement p[ t cf (7) with (3 = 6 dB. 
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